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ABSTRACT 

With a 1999 ASCA observation, PG 2112+059 became notable as the first Broad 
Absorption Line (BAL) quasar found to exhibit a typical radio-quiet quasar X-ray con- 
tinuum underlying a large amount of intrinsic absorption. We present a recent Chan- 
dra ACIS-S3 observation of PG 2112+059 that demonstrates remarkable spectral and 
luminosity variability since that time. In addition to a decrease in the continuum nor- 
malization by a factor of ~ 3.5, the absorption column density has apparently increased 
substantially, and a strong feature in the Fe Ko region has appeared. Concurrent HST 
STIS data compared with archival HST data from earlier epochs show evidence for 
variability of the continuum (up to a factor of ~ 1.7 in the ultraviolet), and in some 
absorption features of the CivA1549 BAL since 1992; however, the Ovi BAL structure 
is consistent with a 1995 observation. We also present evidence for Ly/3-0 VIA1037.62 
and Lya-N VA1242.80 line-locked absorption systems, supporting the assumption that 
ultraviolet line pressure is driving the BAL outflow. Whereas ultraviolet BALs typically 
exhibit only modest equivalent-width variability over timescales of years, the dramatic 
X-ray variability of PG 2112+059 suggests that X-ray spectral variability studies of 
BAL quasars have great potential for probing the physics of quasar winds. 

Subject headings: galaxies: active — quasars: absorption lines — X-rays: galaxies - 
quasars: individual (PG 2112+059) 

1. Introduction 

Recent studies of the strong correlations between supermassive black hole (SMBH) masses 
and the properties of their hosts' galactic bulges (e.g., Ferrarese & Merritt 2000; Gebhardt et al. 
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2000) demonstrate clearly the intimate interrelation between the growth of SMBHs and their host 
galaxies. This connection implies a physical mechanism for regulating the coeval development 
of these structures. During much of their accretion phases, SMBHs reveal themselves as luminous 
quasars (e.g., Yu & Tremaine 2002), and thus quasar winds are likely to provide an important source 
of feedback during the coeval growth of SMBHs and their hosts' galactic bulges (e.g., Fabian 1999). 
These winds are directly observed in the population of Broad Absorption Line (BAL) quasars, which 
exhibits deep, broad absorption lines from high-ionization ultraviolet (UV) resonance transitions. 
BALs arise when our line of sight passes through thick material flowing outward from the nuclear 
region at speeds up to 0.03-0. 3c. Though BAL quasars typically compose ~ 10% of the population 
in optically selected surveys, a recent analysis of the survey criteria for the Large Bright Quasar 
Survey, taking into account selection biases against the BAL quasar population, concluded that 
the true fraction of BAL quasars in the radio-quiet quasar population is 22 ± 4% (Hewett Sz 
Foltz 2003); Reichard et al. (2003) reached a similar conclusion based on the Sloan Digital Sky 
Survey BAL quasar sample. Furthermore, spectropolarimetry, statistics from gravitational lensing, 
and theoretical arguments support this higher intrinsic BAL quasar fraction (e.g., Goodrich 1997; 
Chartas 2000; Krolik & Voit 1998). Therefore, BAL quasars are not exotic objects but rather 
provide an opportunity to study powerful gas outflows that are likely present in typical quasars 
(e.g., Weymann et al. 1991). 

The presence of UV BALs is generally accompanied by extreme weakness in soft X-rays as a 
result of intrinsic absorption with spectroscopically measured column densities typically in the range 
of (1-50) xlO 22 cm~ 2 (e.g., Green et al. 2001; Gallagher et al. 2002). Conversely, X-ray-weak Active 
Galactic Nuclei (AGNs) often show evidence for UV absorption, though not necessarily BALs (e.g., 
Brandt et al. 2000). While significant X-ray and UV absorption invariably occur in concert, the 
absorber properties do not always agree. For example, the column density implied from the X-ray 
observation of the BAL quasar PG 0946+301, N H > 10 24 cm" 2 (Mathur et al. 2000), is significantly 
larger than the best estimates from detailed UV absorption-line studies (-/Vh~ 10 20-22 cm -2 ; Arav 
et al. 2001). Furthermore, among BAL quasars, there are no obvious trends relating the particulars 
of the UV absorption-line properties to the magnitude of the observed weakness in X-rays (Gallagher 
et al. 2003). Such a correspondence between X-ray and UV absorption properties might be expected 
based on theoretical quasar disk- wind models (e.g., Murray et al. 1995) and also on observations of 
Seyfert galaxies where high-resolution X-ray spectroscopy has revealed that the ionized absorbers 
seen in the X-ray regime often share kinematic properties with the UV-absorbing outflows (e.g., 
Kaspi et al. 2002; Gabel et al. 2003). However, at present it is unclear if Seyfert ionized-absorber 
outflows are distinct from or belong to a continuum of winds that extend to the energetic BAL 
outflows in luminous quasars (though see Laor & Brandt 2002 for a discussion of this issue). 

Though a high-resolution X-ray gratings observation of a BAL quasar would offer a direct 
kinematic test of the connection between the velocity structure seen in X-ray and UV absorption, 
the X-ray faintness of this population makes such an observation expensive. Simultaneous X- 
ray/UV variability studies offer an alternative means of investigation. A clear demonstration of 
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absorption variability in one wavelength regime without corresponding changes in the other would 
provide strong evidence that the two absorbers are distinct. Conversely, corresponding changes in 
both X-ray and UV absorption would support an X-ray/UV absorber connection. Understanding 
the relationship between X-ray and UV absorption is essential for constraining important physical 
parameters, in particular the mass-outflow rate from the central engine. The mass-outflow rate 
determines the kinetic influence of quasar activity on star formation in the host galaxy as well as 
the quasar contribution to the enrichment of the intergalactic medium with metals. 

An excellent target for a variability study is the UV-bright, BAL quasar PG 2112+059 which 
has shown significant variability in both the UV continuum and X-ray in the past (Gallagher 
et al. 2001). This quasar is one of the most luminous low-redshift (z = 0.466) Bright Quasar 
Survey objects with My = —26.9 (Hq= 70kms _1 Mpc -1 , Qm = 0.3, and Qa = 0.7 are assumed 
throughout), and it has a broad-band optical polarization of 0.71% (Berriman et al. 1990). HST 
spectra revealed broad, shallow C iv absorption making this a BAL quasar with a balnicity index of 
« 2980 km s^ 1 (Jannuzi et al. 1998; Brandt et al. 2000), well above the lower limit that defines the 
class (Weymann et al. 1991). There are no spectra in the literature that cover the Mgll region, and 
so it is not known if PG 2112+059 has broad Mgll absorption lines (see Boroson & Meyers 1992 
for further discussion of Mgn BAL quasars). However, the lack of Aim absorption (which often 
accompanies Mgll) suggests that PG 2112+059 does not have low-ionization BALs. PG 2112+059 
was detected in soft X-rays with ROSAT (Brandt et al. 2000) which is unusual for a BAL quasar 
(Kopko et al. 1994; Green & Mathur 1996). With the 1999 Oct 30 ASCA data, it was the first 
BAL quasar shown to have an X-ray continuum typical of radio-quiet quasars underlying intrinsic 
absorption. PG 2112+059 showed significant spectral and flux variability between the 1991 Nov 19 
ROSAT and the ASCA observations (Gallagher et al. 2001). 

With the high X-ray flux measured with ASCA as well as its history of X-ray variability, 
PG 2112+059 was a natural target for a coordinated X-ray/UV spectral study to investigate the 
relationship between absorption in the two regimes. In this paper, we present the results from joint 
Chandra/ HST observations of PG 2112+059 taken on 2002 September 1. For reference, the dates 
of the multiwavelength observations presented in this paper are presented in Figure 1. 

2. Observations and Data Analysis 

2.1. Chandra Analysis 

On 2002 September 1, PG 2112+059 was observed with the Chandra X-ray Observatory (Weis- 
skopf et al. 2002) ACIS-S3 detector (Garmire et al. 2003) in 1/8 subarray mode to reduce the 
possibility of photon pile up. The data were processed using the standard Chandra X-ray Center 
pipeline software from which the level 1 events file was used. The data were then reprocessed using 
acis_process_events to remove pixel randomization and filtered on good time intervals, good 
ASCA grades (0, 2, 4, 5, 6), and good status (ignoring the bits indicating afterglow events, which 
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include some source X-rays). The source spectrum was extracted from a circular region with a 4" 
radius; the background spectrum was taken from an encircling annulus with an inner radius of 6" 
and an outer radius of 15". For spectral analysis, the redistribution matrix file (rmf) and auxiliary 
response file (arf) were generated using the standard CIAO 2.2 software, and the arf was then 
corrected for the time-dependent degradation in low-energy quantum efficiency using acis_abs. 6 

Eight keV typically defines the upper limit of the ACIS energy band as the effective area of the 
High Resolution Mirror Assembly drops sharply at this energy; however, in the 7-8 keV band, only 
two counts were within the source region. Therefore PG 2112+059 is not detected above 7 keV. 
We ignore data below 0.5 keV due to calibration uncertainties. We collected 841 photons in the 
0.5-7.0 keV band for a count rate of (1.48 ± 0.05) x 10 -2 ct s _1 . PG 2112+059 was the brightest 
point source detected within the 1'05 x 8^40 ACIS subarray field of view; the next brightest source 
was 2f 27 away with 60+8 full-band counts. The PG 2112+059 radial profile appears consistent with 
that of a point source, and the Chandra position is within 0?5 of the NED 7 optical position. There 
were no background flares, and the quasar showed no significant variability during the 56.868 ks 
observation. 

Based on our best-fitting ASCA model, we had predicted a Chandra ACIS-S3 count rate of 
«9.5xl0" 2 ct s"\ more than a factor of six larger than observed. This large discrepancy is much 
greater than the ASCA/ 'ACIS cross-calibration uncertainties; ASCA and ACIS-S3 fluxes agree 
to within 10% in the 1-10 keV bandpass (e.g., Snowden 2002). Therefore, PG 2112+059 varied 
significantly between the ASCA and Chandra observations. Figure 2a illustrates this point with the 
best-fitting ASCA model convolved with the ACIS-S3 reponse overplotted on the current ACIS-S3 
data. Relative to the ASCA model, the Chandra spectrum appears flatter across the 0.5-7.0 keV 
band pass. Quantitatively, the gross spectral shape can be approximated with a T ~ 0.75 power 
law with the exception of a feature near the FeKa line energy of 6.4 keV evident in the ratio plot. 

As a first attempt to fit the Chandra data, we used the X-ray spectral-fitting software XSPEC 
11.2 (Arnaud 1996) and applied the following model from the 1999 ASCA analysis (Gallagher 
et al. 2001): Galactic absorption (A^h = 6.26 x 10 20 cm -2 ), intrinsic neutral absorption ( Ah = 
[0.86lg;4g] x 10 22 cm -2 ), and a power law (r = 1.95^qH). 8 Even allowing the normalization to 
vary freely did not provide a statistically acceptable fit (x 2 /v = 168.4/47). The intrinsic neutral 
column density as well as the power-law photon index were then allowed to vary. Again, the fit was 
unacceptable (x 2 /v = 84.6/47) with clear systematic residuals near the FeKa region and below 
~ 1 keV in the observed frame. The ASCA data could also be fit with either ionized absorber or 
partial-covering neutral absorber models (though neither was significantly preferred over a simple 

See http:/ /cxc. harvard.edu/cont-soft/software/ACISABS. 1.1. html. 
7 NASA Extragalactic Database, http://nedwww.ipac.caltech.edu. 

8 The ASCA best-fitting model parameters have been adjusted from Gallagher et al. (2001) 
to account for the degradation of the low-energy quantum efficiency of the SIS detectors. See 
http://heasarc.gsfc.nasa.gov/docs/asca/calibration/nhparam.html for further details. 
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neutral absorber), and so we applied the XSPEC models absori and zpcfabs independently, for 
an ionized and partial-covering absorber, respectively. The absorption and continuum parameters 
were free to vary. While each of these fits offered improvements over a simple neutral absorber 
{x 2 / v = 61.6/46 and x 2 l v = 56.5/46 for ionized and partial-covering absorption, respectively), 
clear positive residuals near the FeKa region indicated that another component was required. 
We thus added a cosmologically redshifted gaussian emission line which improved the x 2 value 
significantly for both models. The A% 2 for three additional parameters was 12.5/12.8 for the 
ionized/partial-covering absorber models. In both cases, this indicated a statistically significant 
improvement according to the -F-test at the > 98% confidence level. Given the final fit statistics, 
we slightly prefer the partial-covering absorber model. However, we caution that given the data 
quality, the best-fitting model should be taken as a basic parameterization of the spectrum rather 
than as a truly physical model. In fact, the intrinsic absorber is likely to be much more complex, 
and physical effects such as velocity dispersion and multi-zone ionization structure should be kept 
in mind. The Chandra data, best-fitting model, and residuals are presented in Figure 2b. 

At the 90% confidence level, taking all parameters to be of interest other than absolute nor- 
malization, the best-fitting continuum model parameters are the following: T = 2.09^Qgg and 
jV H = (7.03l|o9) x 10 22 cm -2 , with an absorber covering fraction of / cov = 0. 87to^r The rest- 
frame emission-line parameters are not as well constrained, and so we quote these errors at the 68% 
confidence level: £ linc = 6.38^;^ keV, <7i ine = 0.46^g;^ keV, and EW = 10501)$? eV. Note that 
allowing the width of the emission line to vary versus forcing the line to be unresolved improved 
the fit from x 2 ' l v = 52.3/44 to the final value of x 2 l u = 43.5/43, which is significant at the > 99% 
confidence level. Furthermore, ignoring the observed-frame 4-5 keV region in the spectrum to 
ensure that the emission-line feature was not skewing the continuum fits resulted in best-fit values 
for T, Nn, and / cov that were consistent with those listed above. If this emission line were present 
in the 1999 ASCA spectra with the same flux, it would have been below the detection threshold 
(EW1999 < 383 eV) due to the much higher continuum normalization. The best-fitting spectral 
parameters from the partial-covering absorber model are listed in Table 1; for comparison, the 
values from fitting the same model to the 1999 ASCA data are also presented. The x 2 contours 
of /cov versus iVn for both epochs are shown in Figure 3; though the best-fitting intrinsic column 
density has increased significantly from 1999, f cov remains consistent. This result is upheld with 
simultaneous fitting of the ASCA SIS and Chandra data. 

This intriguing result from the partial-covering absorber model fitting that the intrinsic ab- 
sorption column density had increased significantly over a period of years required additional in- 
vestigation. To that end, we examined whether a simple change in ionization parameter rather 
than column density could account for the spectral changes by simultaneously fitting the ACIS-S3 
and ASCA SISO and SIS1 data with the absori model. We excluded the ASCA GIS data from 
the fits as they are not sensitive below 0.9 keV, the energy regime most important for constrain- 
ing absorber properties. Given that the photon index was consistent between the two epochs, T 
was tied together for all of the subsequent model fitting. The Chandra data between 4-5 keV 
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were ignored to avoid confusion from the Fe feature. For completeness, both the intrinsic column 
density and the ionization parameter, £ = Lx/(nR 2 ) erg cm s _1 , were first tied together, and 
only the continuum normalization was allowed to vary. This resulted in an unacceptable fit, with 
X 2 /v=l5l. 5/117. For Model 1, £ was freed to vary between the Chandra and ASCA epochs with 
a resulting fit statistic of x 2 /v=VSl. 1/116. The residuals for the ASCA data were systematically 
negative below ~ 1 keV, and the best-fitting photon index was unusually flat for a luminous quasar: 
T ~ 1.5. For the subsequent fit (Model 2), the values for £ were tied together for both epochs, and 
the column density was allowed to vary between observations. This resulted in an improvement to 
X 2 /^=123. 2/116 with a more plausible best-fitting T ~ 1.7. Finally, for Model 3, both £ and TVh 
were allowed to vary independently for each epoch with a resulting fit statistic of x 2 /i/=112.0/115. 
The best-fitting parameters from this joint fitting for Models 1, 2, and 3 are presented in Table 2. 
The substantial and significant increase in the intrinsic column density from 1999 to 2002 implied 
by the independent partial-covering absorber fitting is also supported by this ionized absorber joint 
modeling. The change in the low-energy spectrum is not consistent with a pure variation in ion- 
ization parameter. A simple increase in column density is acceptable statistically, though changes 
in both column density and ionization parameter are preferred. Interestingly, with Model 3 the 
Chandra data require a much more highly ionized absorber with £ ~ 500 erg cm s _1 than the 
ASCA data, which are consistent with a neutral gas. Given that the X-ray continuum normaliza- 
tion was significantly higher in 1999, this is contrary to expectation if the X-ray continuum tracks 
the ionizing continuum; the absorber would be expected to be more highly ionized in the presence 
of a stronger ionizing continuum. This suggests that a simple ionized absorber is not appropriate 
in this case. 

The flat spectrum with an effective T ~ 0.75 and strong FeKa emission are reminiscent of 
the reflection-dominated spectra often seen from Seyfert 2 galaxies. However, complex intrinsic 
absorption can fully account for the spectral shape; the best-fitting photon index, T = 2.09+Qgg, 
though loosely constrained, is typical of a direct radio-quiet quasar continuum and consistent with 
the 1999 ASCA data (Figgg = 1.96^o2g)- Furthermore, the latest self-consistent, constant-density, 
ionized reflector models (Ross &: Fabian 1993; Ballantyne et al. 2001) could not provide acceptable 
fits to both the continuum and FeKa emission feature at the same time. Therefore, in the hard 
band we are either directly observing the continuum in a low state or scattered emission from an 
efficient ionized "mirror" that preserves the continuum shape. 



2.2. STIS Analysis 

UV spectra of PG 2112+059 from 1140-3170 A were obtained on 2002 September 1 with 
HST/STIS (Woodgate et al. 1998). Exposure times were 1100 s and 900 s using the G140L and 
G230L gratings, respectively, with the 52" x Qf.'2 slit. The spectra have dispersions of 0.6 A/pix and 
1.58 A/pix, respectively, providing resolving powers of i?=960-1440 for G140L and i?=500-1010 
for G230L across the full wavelength range. The signal-to-noise ratio in the continuum ranged 



- 7- 



from 10 to 25 per pixel. Reductions were made with the standard STIS pipeline, which provides 
flux calibrations accurate to ~ 3% (Brown et al. 2002). The new STIS data are overplotted on the 
previous HST spectra in Figure 4. The 1995 GHRS data covering the O vi emission and absorption 
features appear to be consistent with the current G140L data, while the 1992 FOS data show a 
fainter and redder continuum compared with the current G230L data. 

To investigate potential absorption variability in detail, we focused on the Civ, Lya-Nv, 
and Ovi absorption regions that are not confused with strong emission lines as is the case for 
the Si IV absorption region. First, the data from the previous epochs were convolved with the 
instrumental spread function and resampled to match the lower resolution of the 2002 STIS data. 
Next, wavelength shifts of +1, +2, and +3 pixels were applied to the FOS G270H, FOS G190H, 
and GHRS G140L spectra, respectively, to minimize the differences between previous and current 
epoch spectra. For all but the FOS G190H spectrum, these wavelength shifts were confirmed 
with interstellar absorption lines. To investigate potential variability, the spectra from the 2002 
observation were divided by those from the previous epochs (see the middle panels of Figs. 5-7). 
The division plots illustrate changes in the equivalent widths of absorption features, which should 
appear as significant deviations from the mean level of the ratio plot. The cancellation of broad 
emission lines in the ratio spectra in Figures 5 and 7 indicates that the emission-line changes are 
consistent with their strengths varying in proportion to the changing continuum level. 

In addition, the spectra from the earlier epochs were subtracted from the 2002 spectra as 
shown in the lower panels of Figures 5-7. This comparison is intended to test whether changes 
in absorption features seen in the ratio spectra are caused by the changing continuum level in the 
presence of an unvarying component. For example, BAL troughs frequently contain a contribution 
from an unabsorbed continuum (either from scattering or partial covering of the primary continuum) 
which can fill in the troughs and cause "non-black saturation" (e.g., Arav 1997). In this case, an 
increase in the unabsorbed continuum level while the absorbed component remained constant could 
result in an apparent change in the equivalent width, but such a change would not appear as a 
significant feature in the difference plot. 

In Figure 5 covering the Civ spectral region, the continuum level has increased from 1992 to 
2002 as can be seen clearly in the middle panel. The ratio of /a,stis//a,fos has a median value of 
1.50. In addition, positive deviations from the mean are evident at 1450 and 1481 A at the locations 
of two of the Civ absorption features (marked with asterisks in Fig. 5). The EWs of these two 
features have decreased from 1992 to 2002, at the same time that the continuum normalization has 
increased. These concurrent changes, a decrease in absorption EW accompanied by an increase 
in continuum normalization, are consistent with an increase in ionization that could reduce the 
opacity of individual features. However, these features are not present at a comparable significance 
in the difference spectrum shown in the lowest panel of this figure. This suggests that a constant- 
flux, absorbed component plus a varying, unabsorbed continuum could cause the EW changes in 
the troughs. As mentioned above, this unabsorbed component at the base of the trough could 
arise when an absorber is not physically large enough to cover completely the spatially extended 
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continuum source. The flux from a source other than the direct continuum, e.g., continuum flux 
scattered into the line of sight, could also fill in the absorption troughs. This latter case is perhaps 
more likely, as an unabsorbed component from scattered light would be generated on larger spatial 
scales than the primary continuum and would therefore be expected to respond more slowly to 
changes in the primary continuum. A primary continuum partially covered by the absorber should 
rise and fall in concert with continuum changes. 

The C iv absorption feature seen at rest-frame 1450 A in Figure 5 has the same redshift as the 
1159/1163 A Nv doublet marked with an asterisk in Figure 6. This Nv component shows similar 
behavior: an apparently significant EW change is seen in the ratio spectrum without a correspond- 
ingly significant change in the difference plot. Unfortunately, the Nv absorber corresponding to 
the Civ absorber at 1481 A is not covered by the data in Figure 6. 

Figure 7 with the same plots comparing the O vi spectra from 1995 and 2002 does not demon- 
strate similarly significant deviations in the ratio spectrum. Though the continuum level has de- 
creased by < 25% blueward of the O vi emission line, the absorption structure appears consistent 
between the two epochs, within the uncertainties. If the variability of the C iv absorption EW does 
reflect actual absorber changes, the constant level of the O vi absorption coupled with the changes 
in the Civ structure could indicate absorption variability between 1992 and 1995 with none be- 
tween 1995 and 2002. However, the difference in ionization parameter probed by C iv and O vi 
could also allow a change in C IV absorption with little accompanying change in O vi, depending 
on the ionization structure of the absorbing gas. Unfortunately, the lack of concurrent spectral 
coverage of both C iv and O VI in the earlier epochs does not allow us to distinguish between these 
two possibilities. Furthermore, we have no information on the properties of the UV absorber in 
1999 when the ASCA data were taken. 

In order to investigate the velocity structure in these absorbing features, we also examined 
the GHRS and FOS G190H data at higher resolution. The GHRS data have 0.80 A resolution 
and 0.143 A/pixel, and so the data were binned by 5 pixels to increase the signal-to-noise ratio, 
still giving better resolution than STIS (a 2-pixel resolution of 1.2 A for the STIS G140L grating). 
The FOS data have not been rebinned. The wavelength scale of the GHRS data was corrected by 
+1.00 A, a shift derived from the GHRS wavelength calibration observation. This gave excellent 
agreement with the expected wavelengths of several interstellar lines. 

Attributing the four dips in the broad Civ absorption trough to four C ivAA1548.20, 1550.77 
absorption systems leads to the identification of four Ly/3A1025.72 - O viAA1031.93, 1037.62 systems 
at the same redshifts in the GHRS spectrum. The measured redshifts of the absorption systems 
could then be refined using the higher resolution of the GHRS spectrum, giving z=0.371, 0.387, 
0.402, and 0.418. These are given in Table 3 together with estimated la uncertainties and outflow 
velocities. Figure 8 shows details of the Ly/3 - Ovi, LyaA1215.67 - NvAA1238.82, 1242.80, and 
C iv absorbed regions. We show the positions of absorption lines corresponding to the above four 
redshifts, plotting the data on the same logarithmic wavelength scale so that the four systems show 
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the same relative displacements in the three wavelength regions. The predicted wavelengths for 
Lya and N v are indicated in the middle panel. 

The most striking wavelength coincidences of these four systems, as shown in the middle panel 
of Figure 8, are the overlap of N VA1242.80 in system 1 and Lya in system 3. The same overlap is 
also seen with these lines in systems 2 and 4. System 1 is clearly present in the Lya - N v region, 
while the other systems are less obvious, though all four are evident in O vi absorption in the top 
panel. In the Ly/3 - O vi region, O viA1031.93 of system 1 appears to coincide with Ly(3 of system 
2, and OviA1031.93 of system 3 may also overlap Ly/3 of system 4. 

This type of coincidence, known as absorption line locking, occurs when the relative Doppler 
shifts of gas at different distances from the continuum source are the same as the wavelength 
separation of two strong absorption lines. These systems then become locked into this velocity 
separation, as the system closer to the continuum source absorbs the photons that would continue 
to accelerate the more distant, higher velocity system (e.g., Scargle, Caroff, & Noerdlinger 1970; 
Mushotzky, Solomon, & Strittmatter 1972; Scargle 1973). We suggest here that these outflowing 
systems are locked in relative velocity by the coincidence of N VA1242.80 and Lya lines, and possibly 
OviA1037.62 with Ly/3. Lya-Nv line-locking of sub-troughs of BAL quasars has been discussed 
by Korista et al. (1993) and earlier, for example, by Turnshek et al. (1988) and Weymann et al. 
(1991). If this phenomenon is real it implies that the gas is physically stratified (e.g., system 1 is 
more distant than system 3), material within the BAL region is accelerating, and that radiation 
pressure plays an important role in the dynamics of the outflowing gas (Korista et al. 1993; Arav 
1996). 

The relative absorption strengths of the Ovi doublets in systems 2, 3, and possibly 1, as 
well as the flat bottoms of the Civ troughs which are only partially resolved by FOS, suggest 
that these doublets are close to saturation. Further evidence for saturation may come from the 
Ly/3/Lya absorption ratio in system 1, which appears much closer to unity than to the optically 
thin ratio of 0.16. Given that these absorption troughs are not black, it is probable that there is a 
scattered-light component of the continuum in the O VI region, consistent with the results from the 
absorption-line variability analysis above. Higher resolution and signal-to-noise ratio data could 
well clarify the relationships between the wavelengths and strengths in these UV absorption-line 
systems. Spectropolarimetry would be invaluable for differentiating between the partial covering 
and scattered continuum interpretation. Scattered light is likely to be polarized resulting in some 
polarization of the overall spectrum. In particular, one expects a rise in the polarization level at 
the bottom of the UV absorption troughs as seen in a number of BAL quasars (e.g., Cohen et al. 
1995; Goodrich & Miller 1995; Hines & Wills 1995). This approach holds promise for PG 2112+059 
given its significant low level polarization in very broad band optical polarimetry (Berriman et al. 
1990). 
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3. Results and Discussion 

Between the 1999 ASCA spectra and the 2002 Chandra spectrum presented in this paper, 
the X-ray spectral shape and luminosity of PG 2112+059 have varied dramatically. From the 
best-fitting partial-covering absorber models, the change in the continuum shape at low energies 
can be explained by a significant increase in the intrinsic column density (see Figure 3), while the 
underlying X-ray power-law photon index remained constant within the statistical errors. This 
basic conclusion is not model-dependent; jointly fitting the ASCA SIS and Chandra data with 
ionized-absorber models resulted in the same conclusion (see Table 2). Figure 9 shows the ionized 
and partial-covering absorber models which best fit the data from both epochs. While the power- 
law photon indices are consistent between the two observations, the best-fitting 2002 continuum 
normalization is ~ 30% of the 1999 level. The absorption-corrected, rest-frame 2.0-10.0 keV lumi- 
nosity, -^2. o-io.O) thus dropped by a factor of ~ 3.7 while the absorbing column density increased 
(see Table 1). In addition, structure in the FeKa region is now evident, consistent with a broad, 
neutral emission line with rest-frame EW=580-1570 eV. A feature with this flux would not have 
been detectable in the 1999 ASCA data given the high level of the continuum. 

Comparison of the UV data from 2002 with the two previous epochs, 1992 and 1995, also reveals 
variability, but on much more modest scales than seen in the X-ray. The Civ region observed in 
1992 and 2002 shows evidence for small equivalent width decreases in individual absorption features, 
while the O vi spectra taken in 1995 and 2002 do not show significant evidence for absorption feature 
changes. Over the course of a decade, the UV continuum luminosity varied as well; between 1992 
and 2002 the UV continuum redward of Lya became bluer and brighter. Most of this increase in 
continuum luminosity could have occured between 1992 and 1995 as the 1995 continuum blueward 
of the Lya emission line is within ~ 25% of the 2002 level. 

3.1. Implications of Variability for X-ray and UV Absorbers 

Given the present data quality, translating the significant empirical changes into physical 
understanding is challenging. Unfortunately, there are no UV observations simultaneous with the 
1999 ASCA observations, and so we cannot rule out significant changes in the UV absorption 
properties between 1999 and 2002. However, the concurrent lack of UV absorption-line or large- 
scale continuum variability from the 1995 GHRS to the 2002 STIS observations is worth noting. 
The comparison of these data with the current STIS data is particularly relevant as 1995 is closer 
in time to the 1999 ASCA observation and the continuum level in 1995 is more consistent with 
the 2002 level. As seen in Figure 4, there is no evidence for Ly edge absorption in the 2002 STIS 
G140L spectrum. This puts an upper limit on the column density of neutral hydrogen obscuring 
the UV continuum of Nn < 10 17 cm~ 2 or perhaps a few times larger in the case of partial covering. 
However, the best-fit values for the X-ray absorbing column densities during both the ASCA and 
Chandra epochs are ~5-6 orders of magnitude larger. Either a highly ionized X-ray absorber or 
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an X-ray absorber that does not obstruct the UV continuum could explain this discrepancy. 

Though UV BAL studies are notoriously unreliable for constraining absorber column density 
if partial covering is not properly taken into account (e.g., Arav 1997), a column density increase 
by a factor of ~ 7 as indicated by our X-ray spectral fitting would likely have some effect on the 
UV absorption lines if the X-ray and UV absorption arise in the same gas. The isolated, large-scale 
X-ray variability thus suggests the presence of an X-ray absorber distinct from the UV absorber 
in PG 2112+059. However, without concurrent X-ray and UV data from more than one epoch, 
this claim remains speculative. If the bulk of the X-ray absorption is distinct, it could arise in the 
"shielding gas" postulated in the disk- wind models of Murray et al. (1995). In this picture, the broad 
emission lines and the BALs are produced in the same material at radii of a few light days with the 
shielding gas at smaller distances from the SMBH. This may be typical in BAL quasars; the recent 
detection of relativistically blueshifted X-ray absorption lines in the BAL quasars APM 08279+5255 
and PG 1115+080 (Chartas et al. 2002; Chartas, Brandt & Gallagher 2003), unmatched in velocity 
with UV absorption lines, also supports this picture. 

The four UV absorption systems identified in §2.2 demonstrate, if the line locking is real, that 
radiative forces play a role in the dynamics of the outflow. In particular, Lyman series, O VI, and 
N v line pressure was probably important for accelerating the gas, otherwise the systems would not 
have become locked when these wavelength coincidences occured. However, the Lya-N VA1242.80 
and Ly/3-0 VIA1037.62 line pairs are only two line-locked pairs of many other possibilities among 
UV resonance lines, and so they may only contribute a small part of the total radiation force 
(e.g., Korista et al. 1993). If we speculate that the structure seen in BAL troughs is caused by 
systems locked into a complex network by many absorption-line pairs, then the UV BAL velocity 
structure might be expected to be extremely stable, even in the presence of a varying continuum. 
Confirmation of a change in velocity of an X-ray BAL such as those detected in APM 08279+5255 or 
PG 1115+080 (if line locking plays such a significant role in UV BAL structure) would then suggest 
two intriguing possibilities: (1) the X-ray absorption is occuring in gas as it is being accelerated, 
before the system has become line locked, or (2) the X-ray absorbers are being driven by Compton 
pressure where line locking would not occur. In the latter case, the X-ray absorbers would be 
distinct from the UV BAL gas. 



3.2. Iron Emission 

The combination of the apparent width (a = OA6toll keV) and equivalent width (EW^OSOtf^ eV) 
derived from modeling the Fe emission line in PG 2112+059 is notable. Intrinsically broad lines are 
typically associated with emission from the inner accretion disk; however, the EWs of such features 
are expected to be ~ 150 eV (e.g., Matt et al. 1992). Empirically, the mean for such lines has 
been found to be higher, ~ 220 eV (e.g., George et al. 2000), though still significantly less than we 
observe for PG 2112+059. Therefore, if the line in PG 2112+059 is intrinsically broad, this is hard 
to reconcile with its large EW. Furthermore, a time-delay scenario in which a previously strong 
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continuum has illuminated material whose emission is now more apparent as the continuum level 
has decreased is unlikely with the small scales (~ light hours) of the inner disk. The Fe emission 
from that region would react on short time scales to a diminution of the direct continuum, and so 
a broad, high-EW line from this region would be unlikely. 

In contrast, large EW Fe lines are common features of Seyfert 2 X-ray spectra with Compton- 
thick obscuration to the direct continuum (e.g., George & Fabian 1991; Matt, Brandt, k, Fabian 
1996; Turner et al. 1997). In this case, reflection of direct X-rays off optically thick material usually 
associated with the torus results in a Compton-reflection hump and narrow, fluorescent FeKa 
emission with large EW. Therefore, an alternative picture to an intrinsically broad line is that 
the apparent breadth of the line arises from unresolved, narrow FeKa lines from many ionization 
states. Given the limited statistics and CCD resolution, this scenario is plausible from the positive 
residuals above rest-frame 6.4 keV in Figure 2, though the apparent lack of a strong reflection 
continuum (see §2.1) is problematic unless the optically thick reflector is highly ionized. Given the 
larger spatial scales on which the narrow emission lines are believed to arise, a time-delay scenario 
makes more sense. If this is the case, then the current FeKa emission is reflecting the previously 
strong level of the continuum as was seen in the Seyfert galaxies NGC 2992 and NGC 4051 during 
low continuum states (e.g., Weaver et al. 1996; Guainazzi et al. 1998). In this situation, the FeKa 
line strength would be expected to diminish over time; therefore X-ray monitoring of this source 
has the potential to constrain the spatial scale of any reflecting medium in this luminous quasar. 



3.3. ft x 

X-ray spectroscopy of BAL quasars has indicated that they have typical radio-quiet quasar 
X-ray continua and spectral energy distributions in which soft X-rays are depressed due to intrinsic 
absorption. The power in the X-ray regime relative to the UV continuum is parameterized using 
the quantity a ox =0.384 log(/ !/j 2kcv// 2500A^' wnere ^ ne nux densities, f u , are measured at the rest- 
frame frequencies as labelled. To calculate / 250Q ^> the flux density at rest-frame 2200 A (observed- 
frame 3225 A) was estimated by extrapolating a short distance (~ 55 A) along the continuum from 
the long-wavelength end of the HST STIS spectrum. From that point, the flux density at rest- 
frame 2500 A was calculated using the Sloan Digital Sky Survey composite quasar spectral slope, 
a v = —0.44 (Vanden Berk et al. 2001). We chose this extrapolation to avoid potentially significant 
contamination from a complex of blended Fe n emission lines around 2500 A, and this value for 
the spectral index is consistent with the observed UV continuum. The rest-frame 2 keV (observed- 
frame 1.36 keV) flux density is determined from a local continuum fit to the Chandra data. Based 
on these simultaneous X-ray and UV data of PG 2112+059, the observed value of a ox =— 2.07±0.05 
is within the normal range for BAL quasars (—2.6 to —1.9; e.g., Gallagher et al. 2003; Green et al. 
2001). However, correcting the X-ray spectrum for the intrinsic absorption from the best-fitting 
partial-covering absorber gives a "corrected" a ox ( corr )=— 1.79. In contrast, the value calculated from 
the absorption-corrected 1999 ASCA spectra and the 1992 HST FOS data, a ox(corr ) 1999 = —1.55, 



-13- 



suggested much stronger (a factor of ~ 4.2) X-ray emission relative to the UV continuum, though 
the 1992 UV data were not simultaneous with the ASCA observation. 

As can be seen in Figure 4, the continuum flux densities at rest-frame 2200 A (observed-frame 
3225 A) in both epochs are very similar, and so the decrease in a ox is dominated by the decrease in X- 
ray flux. For comparison, normal radio-quiet quasars of this luminosity typically have a ox ~ — 1.62 + 
0.1 (see eq. 4 of Vignali et al. 2003). The change in a ox ( corr ) of —0.24 over the ~ 3-yr time span is 
thus still within the spread for normal radio-quiet quasars, though a ox (corr) 2002 is at the weak end 
of the distribution. At this point, it is unclear whether the 1999 ASCA spectra or the 2002 Chandra 
data are more representative of the "normal" state of the quasar continuum in PG 2112+059. If the 
1999 ASCA data are more representative, then PG 2112+059 might return to this emission level 
at any time. At least one BAL quasar, PG 1254+047, was argued to be intrinsically X-ray-weak in 
addition to suffering from heavy X-ray absorption (Sabra & Hamann 2001). Though the X-ray data 
quality in that particular case was poor (only ~ 44 photons), it is quite possible that the strength 
of the X-ray continuum of that source varies as well. Probing the timescale and interplay of X-ray 
continuum and absorption variability requires additional monitoring observations of sources with 
demonstrated variability such as PG 2112+059. 



3.4. Comparison of the Bolometric and Eddington Luminosities 

In contrast to the conclusion drawn by Weymann et al. (1991) that all quasars are likely to have 
BAL regions, others have suggested that the BAL phenomenon may instead be a manifestation of 
an evolutionary phase, perhaps of enhanced accretion rate (e.g., Hazard et al. 1984; Becker et al. 
2000). To investigate this issue for PG 2112+059, we estimated the bolometric luminosity as a 
fraction of the Eddington luminosity, Lboi/^Edd, to determine if PG 2112+059 is an outlier in 
the Lboi/^Edd distribution of unabsorbed radio-quiet quasars. From /-band imaging, Fynbo et al. 
(2001) measure / = 18.8 for the apparent magnitude of the luminous elliptical host galaxy which 
corresponds to an absolute magnitude of Mj = —23.2 (in the observed-frame Cousins I band). After 
converting to the rest-frame U-band luminosity (including the appropriate K and color-corrections 
for elliptical galaxies from Table 3 of Fukugita et al. 1995), we obtain an estimate of the black hole 
mass of Mbh=8.8x10 8 M (using eq. 7 from Laor 2001). To estimate the black hole mass from 
the optical spectroscopic properties of the quasar, we utilize -^ 3000 ^ an d H/3 FWHM and eq. 3 in 
Laor (1998) to obtain Mbh=6.8x10 8 M©. Both of these mass estimates have significant systematic 
uncertainties, perhaps as large as a factor of ~ 3 (e.g., Laor 1998; Vestergaard 2002). However, it 
is notable that these two completely independent methods give comparable estimates of the black 
hole mass. 

Using the bolometric correction from the 2500 A monochromatic luminosity for radio-quiet 
quasars of 6.3+3.1 (as calculated empirically by Elvis et al. 1994), we estimate a bolometric lumi- 
nosity for PG 2112+059 of L hoi = (7.7 ± 3.8) x 10 46 erg s" 1 based on L ° measured from the 
2002 STIS spectrum as described in §3.3. Combining these values with the log average of the black 
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hole mass estimates, Mbh=7.7x10 8 M q , gives an Eddington ratio of -^boi/^Edd ~ 0.8. Comparison 
of this rough value with those of similarly luminous quasars in the sample accumulated by Woo & 
Urry (2002, see their Figure 8) indicates that PG 2112+059 does not have an anomalously high 
Eddington ratio at present. Similarly, Yuan & Wills (2003) found that the Eddington ratios of a 
sample of z ~ 2 BAL and non-BAL quasars were consistent. 

3.5. Future Work 

The lack of variability in the velocity structure of UV BALs on timescales of years is one of 
the most serious challenges to the current disk- wind models (e.g., de Kool 1997). While changes 
in the depths of individual absorption-line features have been seen, this EW variability is typically 
ascribed to changing ionization structure within the absorbing gas (e.g., Barlow et al. 1992). Such 
a change in ionization is consistent with the decrease in the depth of C iv absorption features we 
detected during an epoch when the UV continuum increased in normalization as well as became 
bluer. In contrast, changes in velocity of individual absorption components are almost unknown. 
However, the most comprehensive absorption-line variability investigation to date, the 3-yr study 
by Barlow (1993), was not long enough for a conclusive test; the hydrodynamic wind models of 
Proga et al. (2000) generated instabilities on rest-frame timescales of ~ 3 yr. With longer time 
baselines, velocity increases have been claimed for at least one BAL quasar, Q1303+308 (Vilkoviskij 
h Irwin 2001), and an extended study is warranted to investigate this issue further. With high- 
quality UV coverage now extending from 1992 to 2002, almost seven years in the quasar rest-frame, 
PG 2112+059 is an excellent candidate for continued monitoring of its BAL structure. 

The large gaps in coverage of PG 2112+059 over the last decade, illustrated in Figure 1, have 
limited our ability to interpret its multiwavelength variability. Additional epochs of simultaneous 
X-ray/UV monitoring will contribute to understanding the connection between X-ray and UV 
absorption in BAL quasars. The detection of four narrow, probably line-locked absorption systems 
indicates that high resolution UV spectroscopy covering these features would enable a more detailed 
investigation into the partial covering properties of these UV absorption systems, for comparison 
with the partial covering of the X-ray gas. Furthermore, if the bulk of the X-ray absorbing gas 
lies at smaller radii than the UV BAL gas, then the X-ray absorbers might be expected to vary, in 
ionization state, column density, and velocity, on more rapid timescales than the UV BALs. There 
are now examples of BAL quasars with evidence for X-ray absorption variability: PG 2112+059, 
APM 08279+5255, and perhaps PG 1115+080 (Chartas 2000; Chartas et al. 2003). These three 
objects compose a significant fraction of the handful of BAL quasars that have X-ray data adequate 
for variability studies; this demonstrates that the X-ray regime may be a fertile area for exploring 
details of quasar winds inaccessible in the UV. 

Furthermore, to interpret X-ray variability accurately, the hydrodynamic disk-wind models 
(e.g., Proga et al. 2000) may need to be extended to smaller radii (< 10 15 cm) with higher resolution, 
and they must also more explicitly include the effects of Compton and X-ray edge and line pressure. 
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For the highly ionized gas expected at small radii, these will dominate the radiation force over UV 
line pressure (Chelouche & Netzer 2003b, a). More appropriate models for low-energy absorbers 
which include the potentially substantial opacity from velocity-broadened X-ray BALs are also 
required. This line opacity could be comparable to the edge opacity in certain cases. Additional 
items for an X-ray observer's wish list from theoretical modeling are specific predictions for the 
strength, profiles, and variability properties of iron emission and absorption lines. 
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Table 1. Parameters of the Partial-Covering Absorber X-ray Spectral Models 
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iyyy inov 


zuuz sep 


N u (10 22 cm" 2 ) 




o-94±g:52 


' - UO -3.09 


fcov 




n qi+0.09 
u - yi -0.43 


0.87i°; 2 ? 
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1 Q6 +0 - 43 
i.ao_ 29 


2-09±g;S 


keV(obs) (10 -5 phot cnT 


2 s- 1 keV" 1 ) 


29.01 3 ; 2 


8.2l{; 3 






139.8/148 


43.5/43 


Line Energy (keV) 






6.38±°; 42 


a (keV) 






u - 4D -0.32 


EW (eV) 




< 383 c 


1050^° 


^b.5-2.o (10~ 13 erg cm" 2 s" 




2 64+ - 29 


n 26+ - 43 

U - ZD -0.31 


^2.0-7.0 (10~ 13 erg cm -2 s" 




5-84±g;i* 


1 OO+0.18 


£2.0-10.0 (10 44 erg 






1 fiS +0 - 25 
i -°' , -0.17 



a All values listed are for the rest frame, except for fluxes and 
power-law normalizations. Errors quoted are for 90% confidence 
taking all model parameters to be of interest other than absolute 
normalization with the exception of the line parameters, which are 
quoted for 68% confidence. All fits include fixed Galactic absorption, 
A H =6.26xl0 20 cm' 2 . 

h ASCA best-fitting parameters have been adjusted from Gallagher 
et al. (2001) to account for the degradation of the low-energy quan- 
tum efficiency of the SIS detectors (see §2.1). 

c The EW upper limit from the ASCA data is calculated assuming 
the energy and a measured for the Chandra emission line. For an 
unresolved neutral FeKa emission line, the EW< 210 eV (Gallagher 
et al. 2001). 

d Rest-frame 2.0-10.0 keV luminosity corrected for both Galactic 
and intrinsic absorption. 
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Table 2. Ionized- Absorber Model Parameters from Joint Fitting 





ASCA SIS 


Chandra 


Property 21 


1999 Nov 


2002 Sep 



Model 1 

N u (10 22 cm" 2 ) 9.5+£ 2 
£ (erg cm s" 1 ) 4999 b 3414 582+ 533 



r i ™+ - 15 



302 

X 2 /v 131.1/116 
Model 2 

N u (10 22 cm" 2 ) 2.8l 2 ; 3 10.1±f;f 
£ (erg cm s" 1 ) 448±||^ 

1 ± - u '-0.33 

X 2 /^ 123.2/116 

Model 3 

N U (10 22 cm- 2 ) 0.6+ 2 i 11.9+1;| 
£ (erg cm s- 1 ) 0+ 120 502^ 

r i 7 9 +0.25 

1 i -' z -0.24 

X 2 /^ 112.0/115 



a All absorber values listed are for the rest 
frame. Errors quoted are for 90% confidence 
taking 3 parameters (Ah, £, and T) to be 
of interest. All fits include fixed Galactic 
absorption (Ah =6.26x 10 20 cm -2 ), and the 
ASCA SIS spectral models also include the 
corrections for the degradation of the low- 
energy quantum efficiency (see §2.1). 

b This value for £ is "pegged" at the upper 
limit of the allowed range. 
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Table 3. Properties of Line-Locked Absorption Systems 



System 






^outflow 


Number 


^abs 


/^outflow 


( kms -1 ) 


1 


0.37087 ± 0.00017 


0.067 ± 0.004 


20, 100 ± 1300 


2 


0.38725 ± 0.00023 


0.055 ± 0.004 


16, 500 ± 1300 


3 


0.40156 ±0.00010 


0.045 ± 0.004 


13, 500 ± 1300 


4 


0.41837 ± 0.00012 


0.033 ± 0.004 


9, 900 ± 1300 



a The values of /Outflow = v out&ow /c = [(1 + z em ) 2 - (1 + 
z a bs) 2 ]/[(l + Zem) 2 + (1 + <z a bs) 2 ] were calculated using a systemic 
rcdshift of z em =0.466 measured by Boroson & Green (1992) from 
fitting the narrow [O III] emission line. Note, however, that the 
weakness of the [O in] emission in PG 2112+059 introduces some 
uncertainty into this redshift determination, and Corbin & Boro- 
son (1996) measured z cm = 0.460 from broad ~H.fi . The stated 
errors take into account the uncertainties in both z a b s and z em . 
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Fig. 1. — Multi-wavelength light curve of PG 2112+059. Data points are labelled by observatory 
unless otherwise noted. Top panel: X-ray flux densities, in units of 10 -5 phot cm -2 s _1 keV -1 , 
at observed-frame 2 keV from model fits to the local continuum. Middle panel: Dereddened 
UV flux densities, in units of 10 -15 erg cm~ 2 s _1 A -1 , measured at 3225 A in the observed frame 
(2200 A in the rest frame). HST data are indicated by labels with the instrument name. The GHRS 
observation (vertical dotted line) did not cover this region. Bottom panel: Optical V magnitudes 
from the monitoring campaign of Garcia-Rissman et al. (2001). This level of variability is consistent 
with the results from the optical monitoring campaign of Rabbette et al. (1998) who found that 
PG 2112+059 brightened by ~ 0.18 mag in V from 1992 Sep to 1996 Jun. Note that the range of 
the ordinate represents only a factor of two in flux. 
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Fig. 2.— Observed-frame Chandra ACIS-S3 spectra of PG 2112+059. (a) The solid histogram 
indicates the best-fitting ASCA partial-covering absorber model convolved with the ACIS-S3 re- 
sponse. The lower panel shows the ratio of the Chandra data to the ASCA model with a vertical 
dotted line at the energy of the neutral FeKa line; PG 2112+059 has varied significantly in both 
spectral shape and normalization since the 1999 ASCA observation. The Chandra data have been 
truncated at 0.65 keV, the lower limit of the ASCA data, (b) The solid histogram indicates the 
best-fitting partial-covering absorber and broad emission-line model. The ordinate for the lower 
panel, labelled %, shows the fit residuals in terms of standard deviation with error bars of size unity. 
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Fig. 3. — Confidence contours for the two partial-covering parameters, covering fraction (/ cov ) and 
intrinsic neutral column density (-/V H ), in the spectral analyses of PG 2112+059 from the ASCA 
(dotted contours) and Chandra (solid contours) observations. The contours are for 68%, 90%, and 
99% confidence. In both cases, the photon index was allowed to vary freely while calculating the 
contours. The irregularity in the 99% confidence ASCA contour at low Nn is likely due to the 
additional Nn added to the model to approximate the degradation of the low energy ASCA SISO 
and SIS1 detector responses (see §2.1). To determine the confidence contours for the Chandra data, 
the best-fitting values for the emission-line model were fixed. 
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Fig. 4. — Observed-frame UV HST spectra of PG 2112+059. Spectra in black are the current- 
epoch STIS data, while those in color were previously presented in Gallagher et al. (2001). The 
primary broad emission lines are labelled in purple, and the rest-frame Lyman edge wavelength is 
indicated with a solid vertical line. The spectra have been dereddened using the Galactic = 
6.26 x 10 20 cm" 2 to obtain E(B - V) = 0.125. 
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Fig. 5. — Enlargement of the C iv BAL region. The FOS spectrum has been rebinned to match 
the lower resolution of the STIS spectrum. Vertical dotted lines have been drawn at the locations 
of absorption minima to guide the eye. The ordinates for the top and bottom panels are in units 
of 10 -15 erg cm -2 s _1 A -1 . Top panel: The black histogram is the 2002 STIS spectrum, while 
the dashed, green histogram is the 1992 FOS spectrum. Middle panel: Division plot with the 
ratio /a,stis//a,fos plotted at each point. The red histogram is the uncertainty propagated from 
the errors in the individual spectra. Asterisks mark the locations of absorption features discussed 
in §2.2. Bottom panel: Difference plot. The median flux difference is indicated with a horizontal 
line, while the red lines indicate the 1 sigma errors from the median. 
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Fig. 6. — Enlargement of the Lya-Nv BAL region. The FOS spectrum has been rebinned to 
match the lower resolution of the STIS spectrum. Vertical dotted lines have been drawn at the 
locations of absorption minima to guide the eye. Only the wavelength region covered by the G140L 
STIS observation is included; the signal-to-noise ratio of the G270H STIS spectrum is too low in 
this region. Top panel: The black histogram is the 2002 STIS spectrum, and the dashed, green 
histogram is the 1992 FOS spectrum. Middle and Bottom panels: As described in the caption 
to Figure 5. 
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Fig. 7. — Enlargement of the O vi BAL region. The GHRS spectrum has been rebinned to match 
the lower resolution of the STIS spectrum. Vertical dotted lines have been drawn at the locations of 
absorption minima to guide the eye. Top panel: The black histogram is the 2002 STIS spectrum, 
and the dashed, blue histogram is the 1995 GHRS spectrum. Middle and Bottom panels: As 
described in the caption to Figure 5. 
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Fig. 8. — Details of spectra from the three wavelength ranges corresponding (from top to bot- 
tom) to absorption of Ly/3 - Ovi, Lya - Nv, and Civ. The ordinate in each panel is in 
units of 10 -15 erg cm~ 2 s _1 A -1 . The data have been plotted on the same logarithmic scale 
to demonstrate the redshift coincidence of the four systems in these three wavelength regions. 
The SiivAA1393.76, 1402.77 interstellar lines are indicated by dashed lines. There are also several 
unidentified lines, at rest-frame 917.81 A, 920.18 A, 943.17 A, and 952.53 A, the last two being 
indicated by '?'. These unidentified lines could be intervening Lya absorbers. 



- 30 - 




to 

!o L — ■ ■ ■ 

0.5 1 2 5 10 

Observed Energy (keV) 



Fig. 9. — Observed-frame model spectra of PG 2112+059 demonstrating the spectral variability 
between the 1999 ASCA and 2002 Chandra observations. The solid curves are ionized-absorber 
models, and the dashed curves are partial-covering neutral absorber models. All models also include 
Galactic absorption. The emission-line model has the parameters of the best-fitting Chandra model 
(see Table 1); it cannot be positively detected in the ASCA spectra with this normalization. 



